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Most of cosmological observables are light-propagated. I will present coordinates adapted 
to the propagation of null-like signals as observed by a geodesic observer. These “geodesic 
light-cone (GLC) coordinates” are general, adapted to calculations in inhomogeneous 
geometries, and their properties make them useful for a large spectrum of applications, 
from the estimation of the distance-redshift relation, the average on our past light cone, 
the effect of the large-scale structure on the Hubble diagram, to weak lensing calculations. 
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1. Motivations 

The more observational precision increases, the more inhomogeneous our Universe looks. 
SupernovffiP, lensing® and other light-propagated observables will soon encounter asso¬ 
ciated complication^^"^. We present coordinates first developed to simplify averages of 
scalars on our past light-coneM^, next us ed to estimate the effect of inhomogeneities on lu¬ 
minosity distancel^^^^, Hubble diagram^^^l^, and recently applied to lensing quantitie^^ 
(and illustrated in a Lemaitre-Tolman-Bondi model), following this historical order. 

2. The geodesic light-cone coordinates 

We define a light-cone adapted metric (close to “observational coordinates’^^^^, but dif¬ 
ferently composed of 6 arbitrary functions (T, [/“, yab) and totally gauged fixed: 

dsGLC = T^din^ — 2Tdrudr -f 7ab(d6'“ — 17“dw)(d^ — U^dw) . (1) 




Fig. 1. Left: Inhomogeneous light-cone parametrized by GLC coordinates. Right: Homoge¬ 
neous light-cone or an adapted system (like GLC coordinates) in an inhomogeneous geometry. 

This metric uses a null coordinate w defining past light cones, the proper time of a 
geodesic observer r, and angles 0°" that photons keep along their path orthogonal to a 
2-spheres E(ui, r) of constant time in our past light-cone (see Fig. [h . In the FLRW limit: 
w = 7]-f r (conformal time, radius), t = t (cosmic time), {6^,9^) = T = a{t), 

17“ = 0, fab ~ sin^ 6). In general^ T is like an inhomogeneous scale factor (lapse 

function), 17“ like a shift-vector and fab is the metric inside E(w,t). We can notice two 
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direct simplifications in GLC which, combined together, give the distance-redshift relation : 

Redshift: (1 + Zs) = T(wo,ro,?“)/T(wo,rs,0“) , (2) 

Angular distance: = 7 ^/^(sin^^)“^/^ with 7 = det( 7 Qj)) = ^ ■ (3) 

These coordinat es share similarities with historical ones such as “observational coordi- 
nates’^^mSmUn!^ (see elements of comparison in Ref.^^ or the “optical coordinates’^^. 

3. Simplification of light-cone averages 

The light-cone average^of a scalar S (e.g. is in general given by {S) {Vq, Aq) = 

7(5;Vo,^o)/.r(l; Vb, ^ 0 ) and we define the average integral to be gauge invariant, invariant 
under {A, V) — >■ {A{A), V{V)) and general coordinate transformations: 

IiS-,Vo,Ao)= f A‘^xy/^g D{Vo-V)D{A-Ao)M{V,A,d^) S{x) , (4) 

J Mi 


where NiV, A, dfi) is a normalization, D{X) = Sd{X) or 0(A) (Heaviside function): 



Among these 3 types of averages, (>5)^^ is closer to physical observables as it averages 
over the deformed 2-sphere embedded in the light-cone V = Vq and a spatial hypersurface 
A = Aq. In GLC coordinates (where V ^ w, A ^ t) we can simplify the average and use 
Eq. 1^ to get (with Tz = t{zs,Wo, 9°')): 


(5)r 




XiWo,Tz,6^) S{wc 






X{Wo,Tz,ef>) 


(5) 


allowing us to average scalars on the sky, at a certain redshift. 


4. Distance-redshift relation at 0(2) 

The GLC metric enables the computation of d[^{z) at 0(2) in the Newtonian gauge (NG) : 

dsNG = “^(* 7 ) [“(1 + 2'l>)d»7^ + (1 — 24')(dr^ -f r^[A0^ + sin^ 6<d<?!>^])j (6) 

with ^ = Ip -\- , 'i> = Ip + the gauge invariant Bardeen potentials (no matter 

shear at 0(1), or see Ref.^^. Establishing the full transformation between GLC and NG 
coordinates at 0(2) in (scalar) perturbations, we get (t,w, 9^, d'^) = iunc{rj,r,9,(p) and 
(T,17“,7“^) = iunc{tp,ip^'^\ (p^'^'^) which allows us to compute dj^ to 0(2) (using Eq. : 

dLizs,9‘^)^dl^^'^izs)[l + S^s'^izs,9^) + 5‘'^\zs,9^)) . (7) 

The first order in scalar perturbations is given by: 

dW(z„0“) = (^1 - “ Lensing(l) , (8) 
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with T-Ls ~ a'{rjs)/a{ris), Arj = rjo — ris^\ containing (Integrated) Sachs-Wolfe ([I]SW), 
Doppler, and lensing (convergence) effects : 


J = ([a+Q]s - [a+Qlo) - ([9rP]s - [drP]o) ~ SW + ISW + Doppler 
Lensing(i) = = /’'“o, ^ "'"1^ A2'tp{v, Vo - J7,6'“) , 


(9) 


where ‘o’ (‘s’) denote quantities evaluated at the observer (source), and we defined: 


Q{V+,V-,S°')= [ dx ip{ri+,x,e°-) , P{r],r,e°-)= f d??' '^^?/)(r?', r, 6 »“) . (10) 

Similarly obtained 0(2) corrections contain (see Refs.^^^*^ : 

• Two dominant terms : (Doppler)^, (Lensing)^, 

• Combinations of 0(l)-terms: •ips, ([I]SW)^, [I]SWx Doppler, (^,j,Qs)x (Lensing, 
[I]SW, Doppler), 

• Genuine 0(2)-terms : Qs^\ Lensing^^^ = ^ 

• New integrated effects: 4 ^ /’(o)+ da; d+Q + 7 o^ daQ ,x,9°‘), 

• Angle deformations: {'yo)abd+9 °'^^^, 

• Redshift perturbations from Eq. involving transverse peculiar velocities : 

7 ^'’ {daP dbP , daQ dbQ , daQ dbP) , d+ da; [# d+Q + yS^daQ ^ 6 ^] , 
41 da?' ‘^dr [{drPf+yS^daP dbP] , 

• Other important terms : Lens-Lens coupling, corrections to Born approximation. 
These 0(2) results were confirmed recent ly^ ^ work ing directly in terms of GLC coordi¬ 
nates rather than NG. An independent derivation^ ^^ ^ ^^ llead to similar results, but a rigorous 
comparison with Ref. ^ is still lacking. Results were also given for vector/tens or pert ur- 
bation^^ (Poisson gauge) and extended to the case with 0(1) anisotropic stres^^^^^. 


5. Effects of large-scale structure on the Hubble diagram 

In Sec. I^we expressed dj^ in terms of hence it needs a description of the 

Bardeen potentials at 0(1,2). We decompose the first order gravitational potential %p in 
Fourier modes and denote by (...) the ensemble (or stochastic) average over perturbations. 
The 0(2) potentials can be related to ip bylSD: oc v-^{diipd^ip) , diipd^ip. 

Hence we can combine light-cone and stochastic averages in order to study the effect of 
statistical perturbations on the whole sky. For example, the (trivial) average of ip'^ gives: 

¥1) = /o°° '^here P^(fc) = {k^/27v^)\iPk{v)f = A{k/kor^-^T'^{k) {z) 

is the power spectrum describing perturbations. At linear order, with A, Us, fep taken from 
WMAP, T{k) is a transfer function^ including a baryonic component (Silk damping), and 
g{z) is the growth factor describing the recent time evolution of perturbations. In GDM 
we get exactly the spectral coefficients coming from each correction of di^{zs, 9°') described 
in Sec. |^: {di) = ^V^{k) C{kAg). We do the same in AGDM, with reasonable 

assumptions to simplify integration^^, and also with a non-linear power spectrum^^^^^. 
Like d]^ one can also average the flux $ = L/iduSp) ~ 'Fq -f $1 -I- $ 2 - We get: 

• {dl'^) = {df^^^'^)~'^[1 + U{z)] where U{z)zzf{z)f^ V^{k) , 
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• {rfL)( 2 ) = [1 + /d( 2 )] with= -(1/2)/$ + (3/8){($i/$o)^) ■ 

Corrections to involve a flux variance dominated by peculiar velocity and lensing: 

((0i/<^o)")/4 = ((4^')2) ~ (l - { [(([a--P] = )"> + <([a.-P]c)^>] + ((Lensing(i))=^)| . 

as shown in Fig. |^for realistic (non-)linear power spectra^. It turns out that the lumi¬ 
nosity flux is minimally affected by lensing w.r.t. other scalar observables at large redshift. 
This calculation can be seen as a check at 0(2) of Weinberg’s argument of flux conserva¬ 
tion^ and has also been confirmed by recent papers through different approachesl^^^^^Sl 



Fig. 2. Corrections f^{z) and fd(z) for a linear 



(Left) and non-linear (Right) power spectrum. 


Similarly, we can get the average/dispersion of the distance modulus : 

W = 1.25(logio e) [ 2 /$ - ((-I>i/$o)")J , <T^ = 2.5(logioe)\/(($i/cE.o)") . (11) 

Compared to the Union 2 data and using a non-linear power spectrum in ACDM (Fig. 
Left), we And that peculiar velocities explain well the scatter at small 2 : and that lensing 
explains only part of the scatter at large 2 . Finally, we can compare our dispe rsion on 
the Hubble diagram with the experimental estimations coming from lensin^^^®^ (Fig- 
Right). We find that the total effect is well fitted by Doppler (2 < 0.2) -I- lensing (2 > 0.3) 
effects and that the lensing prediction is in great agreement with experiments so far. 




Fig. 3. Left: Theoretical dispersion from power spectrum vs Union 2 data. Right: Dispersion cr^ 
from peculiar velocities and lensing (solid line) compared to experimental estimates (grey areas). 
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6. Jacobi map and weak leasing 

We now consider lensingl^, motivated by Sec. and recent work^^ on galaxy number 
counts in GLC. The relative separation of two neighbour light rays simultaneously emitted 
from a source S and converging to an observer O follows the geodesic deviation eq. (GDE) : 

with = D/dA = the photon momentum, A an affine 

parameter along the photon path, and an orthogonal displacement w.r.t. to the rays. 
We project the GDE on the Sachs basis {s^}a=i ,2 (two zweibeins with flat index Al = 1, 2) : 
3 /u/S^Sb = Sab , ~ 0 , = 0 , H^VxSa = 0 ; with Ufj, the peculiar velocity 

of comoving fluid {S, O comoving too), IlO a “screen” projector orthogonal to and 
M/i + {u°‘koi)~^ kfi. We define the Jacobi map Jg, from the observed sky angle 6^ to = 
by 5"^(A) = Jg{\,Xo)0^. Projected quantities and Rg = Ra0vfik°‘k'^ 
(optical tidal matrix) bring us the Jacobi equation (see e.g. Refs.*^*^ : 

^J^(A,Ao) = R^(A)jg(A,Ao) . (12) 

withl.C.: J^(Ao,Ao)=0 , ^Jii\o,Xo) = {k^^u^,)o5^ . (13) 


A direct resolution of Eq. (12 1 gives the angular distance of the source d^(As) = 
x/dSs = YdetJ^^As,A(^. Also, the (unlensed) angular position of the source 
6^ and the observed lensed position 9^ (of the image) are given by: , 

6^ = /k^u^ , where d-A is the homogeneous and isotropic background our 

model refers to. This allows us to define the so-called amplification matrix as: 


A _ ddf ^ J^(Aa, Ao) _ / 1 - K - 7i -72 -I- w A 

^ ~ ~ dA(Xs) ~V-72-w l-K + jiJ 


which defines the lensing quantities : k — 1 — tr J^/ 2 d^ (convergence), a) = | J 2 — Ji |/ 2 J _4 
(vorticity), I 7 I = = \/(l — k)^ + cu^ — (shear) and fj- = l/(detAl) = 

d^/detjg (magnification). 

Let us now turn to the GLG coordinates and express these lensing quantities in it. 
First, the zweibeins are written as = ( 5 ) 4 , 0 , 3 ( 4 ) and k^ = LuT~^6r (with u a pure 
constant). Second, the solution to Eqs. ( |12[ | and ( |13[ ) is: 

J^(A,Ao) = s^(A) [2u.(7ab)-^]_^sf (Ao) (15) 


where (...) = 9 t(. ..) . The angular distance and the magnification become: 

dA = 2Mro(77o)^^^/A/(det7a{,)o , = $/$ , (16) 


involving d^ = a^(T)r^ with r = w — f dT/a(T) measured from the observer and J? (<I>) the 
flux in the in(hornogeneous) geometry. Expressions for the zweibeins can be obtained in the 
GLC coordinates^, but it is more convenient to compute the squared lensing quantities, 
combined with Sasf = Jab and eAB Sa^b ~ x/l ^ab (e the anti-sym. symbol), to get: 


/(1_^)2+c52A 

V 7? + 7 I ) 


f 


! 

liab~l“'icd 

ad 1 

77 ± 

0 


\dA, 

} 

\ 

(det“'’ 7at,^ 

(det“^7..)J 


(17) 


We thus have general lensing quantities expressed with only 3 metric functions (of 7 af)), 
showing the great advantage of working in GLC coordinates. 
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The Jacobi Eq. ( |12[ ) can be rewritten as a first order differential equation for the 
so-called deformation matrix: 




ai (J2 

<5-2 —CTl 


(18) 


involving the optical scalars, 6 (expansion scalar) and a = ai + (shear scalar), and 
known as the Sachs eqnations : 


^ + |d|2 + 0' = <&oo , ^+2e> = M/o . (19) 

The RHS terms are the Ricci and Weyl focusing and are defined as follows: 

<^>00 = , To = , ( 20 ) 

where Ra.p is the Ricci tensor, Cap^i, the Weyl tensor, and -I- is^- As well as the 

amplification matrix, the deformation matrix simplifies in the GLC coordinates as: 

j A 

(pA a ^ a b ‘ /r»i \ 

= -^Sj^SB^ab ■ (^ 1 ) 

Using s‘\s\ = 7 “^ we get the optical scalars : 


9 = uj 


4T 


oj 7 

4 T 7 


77 _ _ 7 _ det 7 ab 

V4T7y^ 4T2 7 


( 22 ) 


We also get the Ricci and Weyl focusing in GLC coordinates: 

<E>00 = , |To| = ^^(7ac^6d + ^ad^bc _ ^ab^cd-^ y^^Ycd , (23) 

where Tab = 7ab - /^hab - (l/2)7ac7‘'‘^7d6 depends only on ^ab, T and their time 
derivatives. This proves again the usefulness of GLC coordinates for leasing and it was 
illustrated^ by the computation of leasing quantities in the case of an off-center observer 
in a Lemaitre-Tolman-Bondi model (considering only the decaying mode). 


7. Conclusions 

We have shown that there are many advantages in using the GLC coordinates. They are 
indeed adapted to calculations involving light-propagation, they can also be used for weak 
leasing (where 7 ab acts as a screen), and may help to get new predictions on cosmology or 
to study other aspects of leasing (e.g. leasing statistics). 
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